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Objective This study aims to analyze how much heart rate variability (HRV) indices discriminatively respond to age and severity of sleep 
apnea in the obstructive sleep apnea syndrome (OSAS). 

IVIethods 176 male OSAS patients were classified into four groups according to their age and apnea-hypopnea index (AHI). The HRV 
indices were compared via analysis of covariance (ANCOVA). In particular, the partial correlation method was performed to identify the most 
statistically significant HRV indices in the time and frequency domains. Stepwise multiple linear regressions were further executed to exam- 
ine the effects of age, AHI, body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), and sleep parameters on the 
significant HRV indices. 

Results The partial correlation analysis yielded the NN50 count (defined as the number of adjacent R-wave to R-wave intervals differing 
by more than 50 ms) and low frequency/high frequency (LF/HF) ratio to be two most statistically significant HRV indices in both time and 
frequency domains. The two indices showed significant differences between the groups. The NN50 count was affected by age (p<0.001) and 
DBP (p=0.039), while the LF/HF ratio was affected by AHI (p<0.001), the amount of Stage 2 sleep (p=0.005), and age (p=0.021) in the order 
named in the regression analysis. 

Conclusion The NN50 count more sensitively responded to age than to AHI, suggesting that the index is mainly associated with an age- 
related parasympathetic system. On the contrary, the LF/HF ratio responded to AHI more sensitively than to age, suggesting that it is main- 
ly associated with a sympathetic tone likely reflecting the severity of sleep apnea. Psychiatry Investig 2012;9:65-72 
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INTRODUCTION 

Obstructive sleep apnea syndrome (OSAS) is caused by re- 
current obstructions of the upper airway that last for over 10 sec- 
onds during sleep despite breathing effort. OSAS is a common 
disease, with an incidence of 3-7% in adult males and 2-5% in 
adult females.' CUnicaUy, this syndrome has considerable influ- 
ence on personal health, causing excessive sleepiness during day- 
time, narcolepsy, mental disorders including cognitive decline, 
and displeasure or depression, as well as physical complications 
such as obesity, diabetes, gastroesophageal reflux, stroke, pul- 
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monary hypertension, systematic hypertension, coronary artery 
disease, and cardiac arrhythmia." * Impairment of cardiovascu- 
lar system control can also occur in OSAS patients.' Significant 
changes in heart rate variability (HRV) and blood pressure can 
be observed in OSAS patients, which seem to increase the risk 
of cardiovascular disease.*" * 

HRV is physiological variability in the heart rate over time 
due to the interaction between the sympathetic and parasym- 
pathetic nervous systems, and it is influenced by the effect of 
changes in the autonomic nervous system on the sinoatrial node 
caused by both internal and external stimuli.' HRV monitoring 
is a non-invasive method of measuring autonomic nerve func- 
tion in the cardiovascular system. Decreased HRV (i.e., decreased 
variability) indicates abnormal adaptability of the autonomic 
nervous system and is reportedly associated with several dis- 
eases, such as renal failure, hepatic insufficiency, and diabetes, 
as well as cardiovascular conditions such as congestive heart fail- 
ure and myocardial infarction.'" 

Changes in HRV indices were previously suggested as an al- 
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temative diagnostic tool for OSAS." '^ Patients with more severe 
OSAS exhibited a higher low frequency/high frequency (LF/ 
HF) ratio, which reflects the balance between the sympathetic 
and parasympathetic nervous systems.'''" Another study sug- 
gested that the marked rise in very low frequency ( VLF) signal 
during sleep in severe OSAS patients, in contrast to the lack of 
such a change in mild OSAS patients, may represent a tool for 
determining the severity of OSAS.'^ Narkiewicz et al. reported 
that patients with moderate to severe OSAS had shorter RR in- 
tervals (defined as the time between two consecutive R waves in 
the electrocardiogram) and increased sympathetic burst fre- 
quency compared with control subjects."" Furthermore, blood 
pressure variability in patients with moderate to severe OSAS 
was more than double that in controls. Mild OSAS patients ex- 
hibited reduced RR variability compared with controls in the 
absence of any significant difference in absolute RR interval. For 
all patients with OSAS, linear regression analysis showed a pos- 
itive correlation between sleep apnea severity and blood pres- 
sure variability."* However, no dear conclusion can be drawn re- 
garding a correlation between HRV indices and OSAS severity 
due to inconsistent data. 

In addition to OSAS severity, age is an important factor con- 
tributing to HRV changes. In general, regulation of the auto- 
nomic nervous system declines with age, and so HRV indices 
are affected by aging.''''* Vagal activity decreases with age, which 
also affects HRV." This decreased HRV can be explained by the 
fact that the reactions of the autonomic nervous system to exter- 
nal stimulation decrease with age.^° Several HRV indices change 
with age, so that 40% of the SDANN (standard deviation of av- 
erage NN interval) indices and SDNN [standard deviation of the 
NN (normal RR) interval] decrease slowly for about 20 years, 
and 55% of the SDNN index (standard deviation of NN inter- 
val) decreases for about 20 years. Furthermore, pNN50 (the pro- 
portion of NN50 divided by total number of NN) and RMMSD 
(the square root of the mean squared difference of successive 
NNs) decrease more rapidly than all other time-domain HRV in- 
dices.^' One study considering both age and sex suggested that 
both time- and frequency-domain HRV indices associated with 
parasympathetic tone were significantly lower in older than in 
younger men, whereas only the short-term HRV index decreas- 
ed with age in women. This study suggests that the regulatory 
function of the autonomic nervous system as a whole decreases 
in males with age, and so HRV decreases. However, there is no 
significant change in the autonomic nervous system in women 
with age, except in the short-term HRV index. Thus, while HRV 
tends to decrease with both age and increased severity of OSAS, 
changes in HRV indices of OSAS patients reflect the effects of sev- 
eral factors.^^ This study aimed to examine how HRV indices are af- 
fected by OSAS severity and aging in male subjects under the as- 
sumption that HRV decreases with age and more severe OSAS. 



METHODS 
Subjects 

This study was conducted on 176 drug- free adult males who 
were diagnosed with OSAS based on a score of >5 on the apnea- 
hypopnea index (AHI) by nocturnal polysomnography. All 
subjects met the following inclusion criteria: 1) male; 2) >18 
and <60 years old; 3) normal electrocardiogram during wake- 
fidness; and 4) AHI >5. Exclusion criteria were history of intake 
of alcohol or other drugs within 7 days of the nocturnal poly- 
somnographic study, current antihypertensive treatment, pre- 
■vious or current cardiovascular diseases, pulmonary disorders, 
diabetes melMtus, psychiatric disorders, other sleep disorders, 
disorders of the autonomic nervous system or endocrine sys- 
tem that can alter blood pressure, and history of operations or 
continuous positive airway pressure (CPAP) treatment for OSAS. 
Informed consent was obtained prior to the study. Basic infor- 
mation on education level, residence, weight, and stature was 
obtained. Self-reported scales including the Pittsburgh Sleep 
Quality Index (PSQI), Epworth Sleepiness Scale (ESS), Stanford 
Sleepiness Scale (SSS), Morningness-Eveningness Questionnaire 
(MEQ), Beck Depression Inventory (BDI), Beck Anxiety Inven- 
tory (BAI), and State-Trait Anxiety Inventory (STAI) were con- 
ducted to measure symptoms such as insomnia, daytime sleepi- 
ness, depression, and anxiety. All participants provided written 
informed consent, and the protocol was approved by the insti- 
tutional review board. 

Subjects were classified into four groups according to age and 
severity of sleep apnea. Subjects under 40 years old were termed 
young (18-39 years), and those over 40 years old, middle-aged 
(40-60 years). In terms of sleep apnea severity, patients with 
AHI <30 were called mild to moderate (AHI=5-29), and those 
with AHI >30 as severe (AHI=30-120). The young and mild to 
moderate group (n=45), the young and severe group (n=40), 
the middle-aged and mild to moderate group (n=45), and the 
middle-aged and severe group (n=46) were identified as groups 
1, 2, 3 and 4, respectively. 

Polysomnography 

Nocturnal polysomnographic recordings were performed 
using an Embla N7000 system (Medcare-Embla*, Reykjavik, 
Iceland) and Somnologica version 3.3.1 sofhvare (Medcare -Em- 
bla*, Reykjavik, Iceland) during the period when subjects were 
in bed from Ughts-off to lights-on. Electroencephalography was 
monitored using C3/A2 and C4/A1 lead pairs and 01/A2 and 
02/A1 lead pairs. Two pairs of electrooculographic leads were 
used. Electromyographic leads were placed on the submentalis 
and tibialis anterior muscles. Airflow was continuously mea- 
sured by a thermistor and nasal pressure cannula. Respiratory 
movements were monitored using respiratory inductive plethys- 
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mographic belts around the chest and abdomen. Oxygen satu- 
ration was measured by a pulse oximeter sensor on the left sec- 
ond finger 

Sleep stage and events were scored according to the American 
Academy of Sleep Medicine (AASM) manual^' As variables relat- 
ed to sleep structure, time in bed (TIB), sleep period time (SPT), 
total sleep time (TST), sleep latency (SL), sleep efficiency (SE), 
stage 1 sleep (SI), stage 2 sleep (S2), slow-wave sleep (SWS), rapid 
eye movement (REM) sleep, wake, stage 1 and 2 sleep fi-action, 
slow-wave sleep Section, REM sleep fraction, wake fi^iction, and 
arousal index were calculated. The oxygen desaturation event in- 
dex (ODI) was defined as the number of events per hour in which 
oxygen saturation decreased by 4% or more. Hypopnea was de- 
fined as a reduction in airflow by 50-80% for at least 10 s associat- 
ed with an oxygen desaturation of at least 4% or arousal. Apnea 
was defined as an 80% or greater reduction in airflow for at least 
10 s. 24 AHI was calculated by dividing the total number of ap- 
nea and hypopnea events by the number of hours of sleep. 

Data acquisition and analysis 

Electrocardiographic signals acquired by the polysomno- 
graphic machine were digitized at a sampling rate of 250 Hz. 
The RR interval was also used in the analysis. Artifacts were elim- 
inated, and an analysis was performed only for normal beats. 
These processes were conducted using an Embla N7000 system 
(Medcare-Embla* Reykjavik, Iceland) using Somnologica ver- 
sion 3.3.1 (Medcare-Embla', Reykjavik, Iceland). This system 
and software are used in many sleep studies.^"^"^* 

Time-domain variables were mean RR, SDNN, SDNN index. 



RMSSD (defined as the square root of the mean of the sum of 
the squares of differences between adjacent RR intervals), NN50 
count, NN50 of total HR (%), SDANN, and HRV triangular in- 
dex. SDNN is the standard deviation of all RR intervals. RMS- 
SD is the square root of the mean of the sum of the squares of 
differences between adjacent RR intervals. SDANN is the stan- 
dard deviation of the averages of RR intervals in all 5-min seg- 
ments. The SDNN index is the mean of the standard deviation 
of all RR intervals for all 5-min segments. NN50 count means 
the number of pairs of adjacent RR intervals differing by more 
than 50 ms during the entire analysis period. NN50 of total HR 
(%) is the NN50 count divided by the total number of all RR in- 
tervals. The HRV triangular index means the total number of 
RR intervals divided by the maximum height of the histogram, 
excluding boundaries. In the ftequency-domain analysis, the pow- 
er was calculated for very low-frequency (VLF, 0.0033-0.04 
Hz), low-frequency (LP, 0.04-0.15 Hz), and high-frequency 
bands (HP, 0.15-0.40 Hz). The LP/HP ratio was also calculated. 

ANOVA was conducted to compare the demographic and 
polysomnographic data of the four groups. ANCOVA was con- 
ducted to compare HRV indices between the groups, control- 
ling for systolic blood pressure (SEP), diastolic blood pressure 
(DBP), body mass index (BMI), and sleep efficiency (SE). A con- 
trast test was conducted to check differences between groups. 
A partial correlation was performed to select the most statisti- 
cally significant parameter in each of the time- and frequency- 
domain HRV indices, conti-olling for SEP, DBR BMI, and SE. 
Stepwise multiple linear regression analysis was executed for 
the NN50 count and LP/HP ratio that showed the most statis- 



Table 1 . Demographic data of OSAS patients, classified into four groups according to their age and AHI (Mean±SD) 




Group 1 


Group 2 


Group 3 


Group 4 


F 






(N=45) 


(N=40) 


(N=45) 


(N=46) 




P 


Age (years) 


31.29±3.35 


32.75±5.50 


48.51 ±4.69 


49.00±4.75 


194.68 


<0.001 


AHI 


15.62±7.00 


58.38±21.00 


16.42 ±6.66 


54.70±16.70 


121.99 


<0.001 


SBP (mmHg) 


133.31±11.91 


133.90±9.93 


128.47± 14.90 


133.54±11.12 


1.98 


0.118 


DBP (mmHg) 


88.31±10.61 


90.50±12.49 


86.62 ±13.40 


89.93±10.87 


0.94 


0.421 


BMI (Kg/m^) 


25.65±3.26 


28.28±2.79 


25.32±3.50 


27.66±2.85 


9.46 


<0.001 


ESS 


10.11±4.87 


9.82±5.03 


10.58±5.15 


10.93±4.99 


0.40 


0.752 


SSS 


3.23±3.53 


2.94±1.34 


3.45±6.20 


4.88±9.74 


0.63 


0.596 


PSQI 


7.51±3.23 


8.13±2.75 


7.87±3.70 


6.88±3.13 


1.15 


0.331 


BDI 


8.16±6.22 


6.97±4.74 


6.34±5.41 


6.37±5.02 


1.09 


0.354 


BAI 


7.60±6.75 


8.74±7.42 


8.78±8.46 


6.19±6.37 


1.20 


0.312 


STAI-I 


44.75±11.04 


44.55±8.86 


43.81 ±8.22 


45.62±9.40 


0.20 


0.896 


STAI-II 


44.48±9.31 


42.44±8.32 


41.25±9.24 


45.38±9.36 


1.45 


0.231 



Group 1: young, mild to moderate OSAS group, Group2: young, severe OSAS group, Group3: middle-aged, mUd to moderate OSAS group, 
Group 4: middle-aged, severe OSAS group. OSAS: obstructive sleep apnea syndrome, N: number, SD: standard deviation, AHI: Apnea-Hy- 
popnea Index, SBP: systolic blood pressure, DBP: diastolic blood pressure, BMI: Body Mass Index, ESS: Epworth Sleepiness Scale, SSS: Stan- 
ford sleepiness scale, PSQI: Pittsburg Sleep Quality Index, BDI: Beck Depression Inventory, BAI: Beck Anxiety Inventory, STAI: State Trait 
Anxiety Inventory 
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tically significant value in each time- and frequency-domain nificant differences among groups, and of these, NN50 count 

HRV index through partial correlation to examine the effects of showed the greatest statistical significance. Frequency-domain 

age, AHI, BMI, SBP, DBP, and sleep parameter indices (stage 1 
sleep, stage 2 sleep, slow wave sleep, and wake time). Stepwise 
multiple linear regression analysis was conducted with each 
NN50 count and LF/HF ratio as a dependent variable, and age, 
AHI, BMI, SBP, DBP, and sleep parameter indices as independent 
variables. Statistical significance was defined as p<0.05. SPSS 
17.0 was used for data analysis. 

RESULTS 

The mean age of the 176 subjects was 40.65±9.58 years. Table 1 
shows the participants' demographic characteristics. Mean age 
and AHI of each group are shown in Figure 1. BMI (kg/meter^) 
showed significant differences between groups 1 and 2 (p< 
0.002), groups 1 and 4 (p<0.027), groups 2 and 3 (p<0.001), and 
groups 3 and 4 (p<0.006), as determined by the contrast test There 
were no significant differences in PSQI, ESS, SSS, MEQ, BDI, 
BAI, and STAI. 

Variables related to sleep structure and sleep events of each group 
are described in Table 2. Of these, TST, SI, S2, SWS, arousal in- 
dex, ODI, average Sp02, snoring time, and number of limb move- 
ments showed significant differences among groups. 

HRV-related variables for each group are described in Ta- 
ble 3. In the time-domain analysis, RR interval, NN50 count, 
NN50 of total HR (%), and HRV triangular index showed sig- 



Table 2. Polysomnographic data of OSAS patients, classified into four groups according to their age and AHI (Mean±SD) 




Group 1 


Group 2 


Group 3 


Group 4 


F 






(N=45) 


(N=40) 


(N=45) 


(N=46) 




P 


TIB (min) 


408.74±42.10 


413.31147.65 


399.96141.44 


412.67153.48 


0.774 


0.510 


SPT (min) 


397.62±42.30 


398.34146.04 


384.66142.91 


400.01151.52 


1.076 


0.361 


TST (min) 


372.69±47.33 


371.36152.72 


340.88161.90 


359.40162.34 


2.997 


0.032 


SL (min) 


11.12±9.30 


14.97118.33 


15.29118.88 


12.66113.44 


0.728 


0.537 


SE (%) 


92.10±7.I7 


91.0117.48 


86.97113.22 


88.04110.45 


2.635 


0.051 


SI (min) 


56.82±42.02 


78.67154.06 


57.33137.45 


106.38163.05 


9.915 


<0.001 


S2 (min) 


188.56±51.21 


183.04166.46 


171.17142.24 


155.16150.20 


3.527 


0.016 


SWS (min) 


55.10±38.67 


44.76135.37 


35.22128.73 


40.64130.55 


2.841 


0.039 


REM sleep (min) 


71.64±30.25 


62.99130.93 


73.59130.78 


63.22132.00 


1.406 


0.274 


Wake time (min) 


28.84±34.27 


25.23125.75 


41.58143.80 


42.61143.07 


2.367 


0.073 


Arousal Index 


15.44±11.22 


39.18124.15 


14.3019.88 


33.42120.72 


22.708 


<0.001 


ODI 


12.70±6.27 


53.78120.83 


13.7216.61 


50.29118.13 


107.271 


<0.001 


Average SpOz 


95.36±1.50 


92.9012.52 


95.1511.33 


93.0512.75 


17.207 


<0.001 


Snore (min) 


102.42±80.35 


57.38157.94 


76.33169.28 


75.23157.00 


3.288 


0.022 


LM 


32.20±40.03 


91.18194.90 


43.38154.54 


75.59183.19 


6.453 


<0.001 



AHI: apnea-hypopnea index, SD: standard deviation, N: number, OSAS: obstructive sleep apnea syndrome, TIB: time in bed, SPT: sleep peri- 
od time, TST: total sleep time, SL: sleep latency, SE: sleep efficiency, SI: stage 1 sleep, S2: stage 2 sleep, SWS: slow wave sleep, REM: rapid eye 
movement, ODI: oxygen desaturation index, LM: Umb movement 
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analysis showed that VLF power, LF power, HF power, and LF/ 
HF ratio showed significant differences, and of these, HF power 
showed the highest statistical significance, followed by LF/HF 
ratio. 

The mean NN50 count was significantly different among 
groups (p=0.00008). As determined by the NN50 contrast test, 
statistically significant differences were identified between groups 
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Figure 1. Comparison of age and apnea-hypopnea index (AHI) of 
4 OSAS groups. The ages were significant differences between 
young groups (1, 2) and middle-aged groups (3, 4), F=1 94.68, 
p<0.001 , and the AHI were also significantly different between mild 
to moderate OSAS groups (1, 3) and severe OSAS groups (2, 4), 
F=121.99, p<0.001. Group 1: young, mild to moderate OSAS group. 
Group 2: young, severe OSAS group. Group 3: middle-aged, mild to 
moderate OSAS group. Group 4: middle-aged, severe OSAS group. 
OSAS: obstructive sleep apnea syndrome. 
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1 and 3 (p=0.001), groups 1 and 4 (p<0.001), groups 2 and 3 
(p=0.011), and groups 2 and 4 (p=0.001) (Figure 2). RR inter- 
vals were statistically different between groups 1 and 2 (p= 
0.002), groups 2 and 3 (p=0.0002), and groups 3 and 4 (p= 
0.018). The HRV triangular index also showed statistically 
significant differences among groups (p=0.024). Specifically, the 
HRV triangular index was significantly different between groups 

2 and 3 (p=0.007, and groups 2 and 4 (p=0.022), as determined 
by the contrast test. 



Additionally, there were statistically significant differences 
in VLF power, LP power, HF power, and LF/HF ratio, as deter- 
mined by frequency-domain analysis. Specifically, there were 
statistically significant differences between groups 1 and 4 (p= 
0.005) in VLF power, between groups 1 and 2 (p=0.013), groups 
2 and 3 (p=0.004), and groups 3 and 4 (p=0.015) in LF power; 
and between groups 1 and 3 (p=0.001) and groups 1 and 4 (p 
=0.001) in HF power, all as determined by the contrast test. The 
mean LF/HF ratio of was significantly different between groups 



Table 3. Comparison of heart rate variability (HRV) indices between 4 groups in OSAS patients (Mean±SD) 



Domain 


HRV index 


Group 1 


Group 2 


Group 3 


Group 4 


F 


P 






(N=45) 


(N=40) 


(N=45) 


(N=46) 






Time- 


RR interval (ms) 


947. 14± 15.73 


875.20± 16.83 


967.84± 16.03 


913.92±15.38 


5.61 


0.001 


domain 


SDNN (ms) 


108.25±6.38 


109.86±6.83 


88.62±6.50 


96.38±6.24 


2.33 


0.076 


analysis 


SDNN index (ms) 


80.61 ±6.23 


86.07±6.67 


62.82±6.36 


77.46±6.09 


2.27 


0.082 




RMSSD (ms) 


65.36±8.78 


62.14±9.40 


42.70±8.95 


56.64±8.59 


1.23 


0.301 




NN50 count 


6102.24±534.73 


5689.10±572.11 


3597.43±544.97 


3109.26±522.83 


7.68 


0.00008 




NN50 of total HR (%) 


25.10±2.30 


21.58±2.46 


15.63±2.35 


12.64±2.25 


6.00 


0.001 




SDANN (ms) 


66.63±4.27 


59.48±4.57 


57.54±4.35 


56.34±4.17 


1.16 


0.328 




HRV triangular index 


20.21 ±0.91 


21.44±0.98 


17.66±0.93 


18.42±.0.89 


3.22 


0.024 


Frequency- 


Total power (ms^) 


15210.49±905.60 


17166.83±968.90 


15226.81±922.94 


18101.67±885.44 


2.37 


0.072 


domain 


VLF power (ms-) 


8947.23±751.24 


10252.52±803.75 


10075.84±765.63 


11964.33±73453 


2.74 


0.045 


analysis 


LF power (ms^) 


3704.00±278.81 


4745.28±298.29 


3505.48±284.15 


4483.87±272.60 


3.96 


0.009 




HF power (ms^) 


2375.49±175.37 


1990.90±187.63 


1543.81±178.73 


1503.43±171.47 


5.49 


0.001 




LF/HF ratio 


2.13±0.38 


3.20±0.41 


3.01 ±0.39 


4.14±0.37 


4.64 


0.004 



OSAS: obstructive sleep apnea syndrome, SDNN: standard deviation of aU RR intervals, SDNN index: mean of the standard deviation of aU 
RR intervals for all 5-min segments, RMSSD: square root of the mean of the sum of the squares of differences between adjacent RR intervals, 
NN50 count: the number of pairs of adjacent RR intervals differing by more than 50 ms in the entire analysis interval, NN50 of total HR (%): 
the NN50 count divided by the total number of aU RR intervals, SDANN: the standard deviation of the averages of I^l intervals in aU 5-min 
segments, HRV triangular index: total number of RR intervals divided by maximum height of the histogram excluding boundaries, VLF: very 
low frequency, LF: low frequency, HF: high frequency, RR interval: the time between two consecutive R waves in the electrocardiogram, HR: 
Heart rate 
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Figure 2. Comparison of the NN50 count between the four OSAS 
groups. The ISIN50 count was significantly different (F=7.68, 
p=0. 00008) when adjusted by systolic blood pressure (SBP), dia- 
stolic blood pressure (DBP), body mass index (BMI) and sleep effi- 
ciency (SE). The contrast test further showed that Group 1 was 
significantly different from Group 3 (p=0.001) and Group 4 
(p<0.001). Similarly, Group 2 was significantly different from Group 
3 (p=0.011) and Group 4 (p=0.001). OSAS: obstructive sleep ap- 
nea syndrome, NN50 count: the number of pairs of adjacent RR 
intervals differing by more than 50 ms in the entire analysis interval. 



Figure 3. Comparison of the LF/HF ratio between the four OSAS 
groups. The ratio was significantly different (F=4.64, p=0.004) 
when adjusted by systolic blood pressure (SBP), diastolic blood 
pressure (DBP), body mass index (BMI) and sleep efficiency (SE). 
The contrast test further showed that Group 1 was significantly dif- 
ferent from Group 4 (p<0.001). Similarly, Group 3 was significantly 
different from Group 4 (p=0.040). OSAS: obstructive sleep apnea 
syndrome, LF/HF ratio: Low frequency/High frequency ratio. 
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Table 4. Partial correlation of age and apnea-hypopnea index (AHI) with heart rate variability (HRV) indices controlling systolic, diastolic 
blood pressure, body mass index and sleep efficiency of OSAS patients 





Variable 




Age 






AHI 




ip 


P 


Df 


ip 


P 


Df 


Time-domain analysis 


RR Interval (ms) 


0.084 


0.2717405 


170 


-0.304 


0.0000518 


170 




SDNN (ms) 


-0.220 


0.0037716 


170 


0.013 


0.8687622 


170 




SUNN mdex (ms) 


-0.164 


0.0320747 


170 


0.093 


0.2239395 


170 




RMSSD (ms) 


-0.108 


0.1581157 


170 


0.005 


0.9462984 


170 




NN50 count 


-0.373 


0.0000005 


170 


-0.058 


0.4502996 


170 




NN50 of total HR (%) 


-0.337 


0.0000062 


170 


-0.098 


0.1987834 


170 




SDANN (ms) 


-0.194 


0.0106428 


170 


-0.121 


0.1127152 


170 




HRV Triangular index 


-0.319 


0.0000198 


170 


0.072 


0.3467232 


170 


Frequency-domain 


Total power (ms^) 


-0.001 


0.9925901 


170 


0.260 


0.0005751 


170 


analysis 


VLF power (ms^) 


0.107 


0.1621830 


170 


0.212 


0.0051794 


170 




LF power (ms^) 


-0.076 


0.3205831 


170 


0.356 


0.0000016 


170 




HF power (ms^) 


-0.309 


0.0000378 


170 


-0.147 


0.0547991 


170 




LF/HF ratio 


0.230 


0.0023515 


170 


0.432 


0.0000000034 


170 



OSAS: obstructive sleep apnea syndrome, RR interval: the time between two consecutive R waves in the electrocardiogram, SDNN: standard 
deviation of all RR intervals, SDNN index: mean of the standard deviation of all RR intervals for all 5-min segments, RMSSD: square root of 
the mean of the sum of the squares of differences between adjacent RR intervals, NN50 count: the number of pairs of adjacent RR intervals 
differing by more than 50 ms in the entire analysis interval, NN50 of total HR (%): the NN50 count divided by the total number of all RR in- 
tervals, SDANN: the standard deviation of the averages of RR intervals in all 5-min segments, HRV triangular index: total number of RR in- 
tervals divided by maximum height of the histogram excluding boundaries, VLF: very low frequency, LF: low frequency, HF: high frequency 



Table 5. Stepwise multiple linear regression analysis of variables' 
relationships to the NN50 count and LF/HF ratio in OSAS patients 



Step Variable 


R 


Adjusted R^ ^ 
R^ change 


SB 


t P 






NN50 count 






1 Age 


0.38 


0.14 0.14 -149.47 


27.13 


-5.51 <0.001 


2 DBF 


0.41 


0.16 0.02 -45.53 


21.91 


-2.08 0.039 






LF/HF ratio 






1 AHI 


0.46 


0.21 0.21 0.05 


0.01 


6.65 <0.001 


2 S2 


0.51 


0.25 0.04 -0.01 


0.003 


-2.87 0.005 


3 Age 


0.53 


0.27 0.02 0.04 


0.02 


2.33 0.021 



NN50 count: the number of pairs of adjacent RR intervals differing 
by more than 50 ms in the entire analysis interval, LF/HF ratio: Low 
frequency/high frequency ratio, OSAS: obstructive sleep apnea syn- 
drome, SB: standard errors, DBF: diastolic blood pressure, AHI: Ap- 
nea-hyponea index, S2: stage 2 sleep 

1 and 4 (p<0.001) and groups 3 and 4 (p=0.040), as determined 
by the contrast test (Figure 3). 

By partial correlation, controlling for SBP, DBP, BMI, and 
sleep efficiency, the parameters with the greatest statistical sig- 
nificance for each time- and frequency-domain HRV index were 
identified as NN50 count (rp=-0.373, p=0.0000005) and LF/HF 
ratio (rp=0.432, p=0.0000000034), respectively (Table 4). To 
determine which variables had the greatest effect on NN50 
count and LF/HF ratio, a stepwise multiple linear regression 
analysis was conducted (Table 5). Of age, AHI, BMI, SBP, DBP 



20,000 - 



15,000 - 



R^=0.14 



K 10,000- 



5,000 




10.0 20.0 30.0 40.0 50.0 60.0 
Age 

Figure 4. The linear regression between age and NN50 count of 
the four OSAS groups. The NN50 count was affected by age (p< 
0.001) and diastolic blood pressure (p=0.039) in the order named 
in the stepwise multiple linear regression analysis. LF: low fre- 
quency, HF: high frequency, NN50 count: the number of parts of 
adjacent RR intervals differing by more than 50 ms in the entire 
analysis interval, OSAS: obstructive sleep apnea syndrome. 

SI, S2, SWS, and wake time, the NN50 count was most affect- 
ed by age (p<0.001) and DBP (p=0.039), and the LF/HF ratio 
was most affected by AHI (p<0.001), stage 2 sleep (p=0.005), 
and age (p=0.021), in the order named (Figures 4 and 5). 
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Figure 5. The linear regression between AHI and LF/HF ratio of the 
four OSAS groups. The LF/HF ratio was affected by AHI {p<0.001), 
Stage 2 sleep (p=0.005) and age {p=0.021) in the order named in 
the stepwise multiple linear regression analysis. LF/HF raito: Low 
frequency/high frequency index, AHI: Apnea-Hyponea Index, OSAS: 
obstructive sleep apnea syndrome. 

DISCUSSION 

Of the time-domain HRV indices, RR interval, NN50 count, 
NN50 of total heart rate, and HRV triangular index showed sta- 
tistically significant differences, and of the frequency-domain 
indices, VLF power, LF power, HP power, and LF/HF ratio show- 
ed statistically significant differences among the four groups, 
according to age and the severity of OSAS. NN50 count and LF/ 
HF ratio tended to decrease with age and with increased AHL 
This result confirms those of previous studies." "'^''^' Our data 
suggested that NN50 and LF/HF ratio responded differently to 
aging and AHI severity. NN50 count decreased more with age 
than with AHI severity, and the LF/HF ratio increased more with 
AHI severity than with age. 

Thus NN50 count and LF/HF ratio showed the greatest sta- 
tistical significance. Therefore, we determined the effect of sev- 
eral variables on NN50 count and LF/HF ratio by stepwise mul- 
tiple linear regression analysis. Only age had a significant effect 
on NN50 count. In the case of LF/HF ratio, AHI was found to 
have the greatest effect, followed by S2 and age. 

The NN50 count indicates how often the RR interval differs 
by over 50 ms from the average RR interval, and it strongly re- 
flects the activities of the parasympathetic nervous system. In 
general, the fionction of the parasympathetic nervous system is 
known to decrease with age.™ Therefore, NN50 count will de- 
crease with age. This is consistent with our data; age showed the 
greatest correlation with NN50 count, and only age was found 
to have significant affect in the stepwise multiple Unear regres- 
sion analysis. 

LF/HF ratio reflects the overall balance of the autonomic 



nervous system. In some cases, this ratio is used as an index of 
activity of the sympathetic nervous system. Hypercapnia and 
hypoxia occur in OSAS patients due to repeated sleep apnea, 
which, in turn, increases sympathetic activity through other me- 
chanisms, including the chemoreceptor reflex, and accentuates 
sympathetic nerve activity.''''^ Sympathetic nervous system ac- 
tivity is known to increase with age.'^ In this study, the LF/HF 
ratio increased with AHI and age, which is consistent with the 
fact that sympathetic nervous system activity increases with both 
age and AHI and that the LF/HF ratio reflects the sympathetic 
nerve activity. The LF/HF ratio showed the highest correlation 
with AHI, and in the stepwise multiple linear regression analy- 
sis, the LF/HF ratio was found to be most affected by AHI. 
Therefore, the LF/HF ratio, which is closely related to sympa- 
thetic nervous system activity, may be affected more by the se- 
verity of sleep apnea than by aging, which confirms the conclu- 
sions of previous studies." S2 had a greater effect on the LF/ 
HF ratio than did age, perhaps because S2 increases with AHI 
severity^'' 

Although this study used a large population (n=176), it has 
several limitations. First, female subjects were not included, and 
so changes in female subjects could not be examined. Howev- 
er, this can also be regarded as an advantage, as the study pop- 
ulation was more homogenous. In a previous study, time- and 
frequency-domain analyses associated with aging showed many 
significant changes in male subjects, but only the short-term 
HRV index changed in female subjects.^^ Therefore, inclusion 
of only male subjects may lead to a superior evaluation of chang- 
es in HRV. Second, the first-night effect was not considered in 
this study. Increased sleep latency, decreased overall sleep time, 
increased wake frequency, and invasion of alpha in non-REM 
sleep may occur due to first-night effects.'^^ However, polysom- 
nography was conducted under the same conditions on all sub- 
jects, and thus first-night effects likely did not influence the 
comparisons among groups. Third, this study did not compare 
OSAS patients with normal controls. Thus, the ability to deter- 
mine whether changes in HRV indices were quantitatively mean- 
ingful compared with controls was limited. 

This study analyzed the effects of age and OSAS severity on 
HRV indices and indicated that these effects differed among 
HRV indices. The LF/HF ratio, which mainly is associated with 
sympathetic tone, is more likely to reflect changes in sleep ap- 
nea severity. NN50 count, which reflects mainly the parasym- 
pathetic system, is more likely to reflect changes in age. There- 
fore, our data suggest that age and severity of sleep apnea exhibit 
different effects on various HRV indices. 
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